Development of an approach to assess critical 
doses and dose rates for cultivated plants 


A.A. Oudalova A., L.N. Ulyanenko and S.A. Geras’kin 


Russian Institute of Agricultural Radiology and Agroecology, 249032 Obninsk, 
Russia 


Abstract. The aim of this work is to develop methods for an assessment of 
critical doses and dose rates that can result in significant radiation-induced 
effects in agrocenoses. This is realized on an example of cultivated plants which 
are one of the main components of agroecosystems. Available information on 
dose dependences in such umbrella endpoints as reproduction, survival, 
morbidity, alterations in morphological and biochemical processes, genetic 
effects in crops, vegetables, fruit trees, etc are gathered from papers issued 
mainly in Russian scientific press during last 50 years. Data are maintained as 
database in MS Access that contains about 10000 records at the moment. As 
critical ecotoxicity values, doses are considered producing 50% changes of 
biological effect at acute impact (EDso), or dose rates resulting in 10% changes at 
chronic exposure of plants (EDRio). The critical doses and dose rates for 
agrocenosis estimated from available information on reproduction and survival 
are presented here. Primary data are assesed for their quality according to 
several criteria. Three models (linear, logariphmic, logistic) for dose-effect 
relationship are tested for an applicability to fit a dose-effect dependence taking 
account of their goodness-of-fit and robustness of EDs) and EDRio estimates. 


1. INTRODUCTION 


In recent years many efforts have been undertaken to develop a system of 
radiation protection for non-human biota [1, 2]. Agrarian ecosystems are of 
special concern from the viewpoint of establishing safe levels of radiation 
impact on the environment because a) their contamination can affect human 
health via radionuclide uptake with food, and b) agroecosystems are ones of the 
most sensitive biocenoses to a number of environmental impacts including 
ionising radiation [3]. While the existing system of radiation standards [4] is 
effective in providing radiation protection of human and restricts radionuclides 
content in food chains, there are no guidelines on setting any limitations to 
directly protect agrarian ecosystems from negative effect of radiation or to 
demonstrate their protection at radiation levels safe for human. The aim of this 
work is to develop methods for an assessment of critical doses and dose rates 
which violation can result in significant radiation-induced effects in 
agrocenoses. This is realized on an example of cultivated plants as the main 
component of agroecosystems. Cultured plants not only contribute essentially to 
food production but also fulfil an important ecological function in agrarian 
biocenosis. 


2. MATERIALS AND METHODS 


2.1 Database 


Available information on radiation-induced effects in cultured plants are being 
gathered from papers issued mainly in Russian scientific press during last 50 
years; the work is ongoing. Data are maintained as database in MS Access that 
now contains about 3000 records from 175 original sources (161 - in Russian). 
There are available about 10000 datasets which are pairs of “exposure level - 
biological effect” data accompanied with information on experimental design, 
species and its life stage, exposure parameters and conditions. Quantitative 
data are collected for about 60 species of cultivated plants. Most information 
found concern radiation effects in crops, vegetables and roots, as well as 
legumes - 38%, 21% and 18% of all entries, correspondingly. 

Biological responses are grouped under umbrella endpoints of 
reproduction, survival, morbidity, alterations in morphological and biochemical 
processes, and (cyto)genetic effects. Only analysis of reproduction and survival 
endpoints at acute and chronic irradiation of plants during vegetation are 
presented in this paper. Reproduction endpoint pools traits of biological 
productivity (seed and straw mass, dry and wet biomass, number of plants per 1 
m’, etc) and reproductive potential (portion of fertile pollen seeds, 
fertility/sterility, number or portion of male and female flowers, etc). The 
survival indexes are germination ability and rate (in field or in laboratory), 
survival of germinated seeds or plants by certain time/day, etc. 


2.2 Dose (rate) - effect models 


Each dataset was substituted for regression analysis to reconstruct the 
dose(rate) - effect relationship and then estimate the critical ecotoxicity values 
that were defined the same way as in [5, 6]. Thus, the dose giving 50% change 
in observed effect (EDs.) in comparison to control is considered as the critical 
dose for acute exposure, and the critical dose rate for chronic exposure is 
defined as the dose rate resulting in 10% change of observed effect (EDRio). 
Three models were tried and tested for an applicability to fit the dose(rate)- 
effect dependences taking account of their goodness-of-fit and robustness of 
EDs 9 and EDRio estimates: 


A. linear model y=at bx, 
B. logarithmic model y=a+ bln(x, 
x 
C. logistic model =q 1- —— |. 
j 4 riy 


For the logarithmic model spontaneous effect data at x = 0 can not be used 
because the logarithm argument must be over zero (x > 0). 


2.3 Data quality assessment 
At the EDs59ÆDRıo computation, primary data were assesed for their quality 


according to several rules, some of them were similar to those proposed in [5, 6] 
for data quality assurance: 


Rule 1. Volume of dataset n must be enough for regression analysis with a given 
model, i.e. degrees of freedom dfshould be >0, where df= n-np-1, np - the 
number of free parameters in a model. 

Rule 2. The variation of effect with dose or dose rate should be monotonous and 
consistent with the current knowledge in radiobiology. 

Rule 3. Regression should be significant at least at p < 0.10. 

Rule 4. At least one effect observation Z should be located within 10 to 90% of 


the effect variation range [Yx=0, Yx= ©]; 

Rule 5.The estimated critical ecotoxicity value should be within the dose(rate) 
range studied: 
Xmin < EDs5o < Xmax OF Xmin < EDRio < Xmax. 


3. RESULTS AND DISCUSSION 
3.1 Datasets of good quality 


In total there were available 498 datasets on the endpoints of reproduction and 
survival. Over 80% of them did not meet the requirements of the five rules to 
the data quality (Table 1). In general, the number of proper datasets depends on 
the dose-effect relationship chosen because of different degrees of freedom and 
varying goodness-of-data fitting. At the application of the linear and logistic 
models less data were rejected then with the logarithmic model. Many original 
studies give only two observations typically referring to the control and 
impacted groups. This is why up to 70% of the datasets are rejected by Rule 1. 
Among the rest data, many datasets do not meet requirements of sufficient 
quality of data fitting, i.e. regression significance (Rule 3). Also, Rule 5 
appeared to be quite restricting for acute exposure data, and Rules 2 - for 
chronic exposure (likely because of observations of stimulating effect at low 
dose rates). 


Table 1. Quality assessment for data on reproduction and survival of cultured plants at 
acute and chronic exposure during vegetation 


Chronic exposure (168 
datasets) 
Number of datasets meeting requirements 


Acute exposure (330 datasets) 


Rule 


Model 


Linear (np = 2) 


Logarithmic (np 
) 


Logistic (np = 3) 


3.2 Critical ecotoxicity values 


Fig. 1 illustrates an approximation of data on reproduction endpoint by three 
dose(rate)-effect models and estimates of the critical dose and dose rate at 
acute and chronic exposure, correspondingly. Calculated critical ecotoxicity 
values are presented in Table 2 (acute exposure) and Table 3 (chronic exposure) 


as average EDso (EDRio) with standard error and range of variation for a given 
species. 

Data of good quality on survival are scarce both at acute and chronic 
exposure. Critical ecotoxicity values are obtained only from one or two datasets 
for every species. For few species, where data are available for both endpoints 
considered (corn, rice, soya for acute, wheat and barley for chronic exposure), 
EDso and EDRio are lower when estimated from reproduction (Table 2 and 3). 
Thus, reproduction seems to be more radiosensitive criteria than survival, which 
is in line with our previous findings [9] and other studies [11, 12]. 
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Figure 1. Approximation of data with three models and estimation of the critical dose (A) 

and dose rate (B) 
A - Yield of grains from corn plants exposed on the 54 day of growing [7]; EDso =15.9, 

15.8, and 15.5 Gy from linear, logarithmic and logistic models, correspondingly. 
B - Fertility of the main ear of wheat plants exposed in y-field from 2-week seedlings to 

harvesting 16 h per day [8]; EDRio = 5.40, 7.92, and 10.1 mGy/h 

It is obvious that critical ecotoxicity values depend on the model used. An 
analysis of the data available shows that at acute exposure the differences in 
EDs 9 obtained from three models appeared to be insignificant (p > 0.05) for all 
species both for reproduction and survival endpoints. 

In the case of chronic irradiation, the critical dose rates EDRıo obtained with 
the logistic model are always higher than that from the linear one (Table 3), 
sometimes up to 2 orders of magnitude. This is because a logistic curve changes 
very slowly at low dose rates, and thus reaches the 10% decrease from the 
control level at higher values of dose rate (Fig. 1B). But variance of the ‘logistic’ 
EDRip values estimated from different datasets for the same species is so high 
that the average EDRio does not differ from 0 for most cases (Table 3). This 
makes differences with the ‘linear’ EDR; values insignificant in all cases. Such a 
high variability is bound to a lack of good quality effect data, which is 
recognized for chronic radiation exposure of biota in general [13] and 
agricultural plants, in particular [9]. Also, the original studies differ very much 
in methods and study design, dosemetric information is often incomplete, 
studied ranges of dose rates poorly overlap. 


Table 2. Critical doses for species of cultivated plants after acute exposure 


Linear Logarithm Logistic 
N EDso, Gy N EDso, Gy N EDso, Gy 


Species 


mean+s | min-max mean+ min- mean+s min- 


e | | se max e max 


Reproduction 

Corn 5 11.91. ; 13.42. 
2 

Oats 5 20.945. ý 15.62. 
4 

Wheat RR roe i 9.9+1.4 

Rice 2 219.3+ 211-228 118+19 

Barley 6| 9.9+1.2| 5.3141 

soya a E 

Potato 7 54.11 22.5-154 

Strawberry | 1 101.5 

Clover 3 222+64 137-348 

Fescue 2 33.6+4 29.6-37.7 

Ryegrass 2 21.342. 19.2-23.5 

Cotton 1 66.7 

Corn 

Rice 1 244 

Beans 2 | 1.9+0.4 1.48-2.33 

Soya 2 | 32.8+1. 31.1-34.5 

Agropyron 1 21.8 

Tobacco 2 | 0.02+0. 

plant 


N - number of datasets 


3.3 Comparison to previous findings 


In our previous contribution [9], the critical doses to species of cultivated plants 
were estimated from the linear dose-effect relationship. Original effect data 
were transformed in relative units as % of corresponding controls, pooled 
together to form arrays under every specific endpoint, and submitted for 
regression analysis. That way of data handling does not require rejecting 
datasets that violate the rules of data quality and, thus, allows using much more 
data for calculations than the approach used in present study. Indeed, the 
pooled arrays included together nearly 600 pairs of data on reproduction and 
survival [9], while 81 datasets meeting the five rules requirements for linear 
model in this work (Table 1) contain only about 370 paired data entries even 
despite the database has been enlarged from 7000 to 10000 entries. 

Findings from two approaches (pooled arrays and individual datasets) are 
compared in Table 4. Critical doses or doses rates calculated from reproduction 
data are presented when there are estimates for the same plants. To ensure 
some confidence in critical doses estimated from pooled data, two requirements 
were fulfilled in [9]: 1) regression significance of at least 10%, 2) a minimal 
number of datasets N = 10 for single species. In present study it is considered 


that averaged values for EDs) and EDRio are not safe enough if they are based 
on less than 3 datasets (N < 3). From Table 4, if the linear approximation is 
applied, there are no significant differences between critical ecotoxicity values 
obtained from pooled and individual datasets for all species except ED; for 
barley. The approximation of data after acute exposure with the logistic function 
gives typically the lower estimates than in the first approach, and this difference 
is significant for oats and wheat. However, in the case of chronic exposure 
variance of EDRio calculated from logistic dependence is extremely high so that 
EDR, does not differ from 0. At the same time, linear extrapolation gives 
reasonable values from both approaches. Their comparison is possible only for 
one species - wheat, and the critical dose rates from pooled and individual 
datasets are nearly equal. 


Table 3. Critical dose rates for species of cultivated plants after chronic exposure 


Species Linear Logarithm Logistic 
p N EDRio, mGy/h N EDRio, mGy/h N EDRio, mGy/h 
mean+ min- mean+ min- 
se max se max mean+se | min-max 
Reproduction 
Wheat 5.4- 9.7+1. 7.9- 9.8- 
7 | 6.5+0.8 11.2 2 5 11.0 9 | 230+218 1977 
Barley 139412 126- 12058. 
2 9 152 1 3.8 3 2 4.0-185 
Bean 3| 8.01.7] j74 2| 7*0 | 6173] 3| 16.026.1 | 9.0-28.1 
Peas 523+43 38- 1020+77 129- 
3 ol 1380| 1| 1018 E 5 2567 
Currant 1 29.2 1 5.5 
Lettuce 2030+20 27.8- 
1 18.5 1 4.7 2 00 4030 
Tomato 
plant 1 6120 
Tobacco 1 4.7 1 3.7 1 8.5 
plant 
Survival 
Wheat 1 3640 1 3660 
Barley 1 2150 
Lettuce 1 5750 1 5760 
Tomato 1 7030 1 7090 
plant 
Clover 1 12740 1 12800 


Table 4. Critical ecotoxicity values with 95% confidence intervals within brackets 
obtained from pooled [9] and individual (this paper) datasets for reproduction 


endpoint 
Dataset Pooled Individual 
Model linear linear logistic 
ED5o, Gy (acute exposure) 
Oats 27.8 (24.4+32.5) 20.9 (5.6+36.1) 15.6 (9.4+21.7) 
Wheat 18.3 (16.4+20.4) 12.5 (9.2+15.7) 9.9 (6.7+13.1) 
Barley 15.6 (14.6+16.6) 9.9 (6.7+13.1) 9.9 (4.0+15.7) 


Potato 45.8 (37.8+42.7) 54.1 (8.6+99.6) 39.6 (0+94.5) 
EDRio, mGy/h (chronic exposure) 


Wheat 6.8 (3.3+10.0) 6.5 (4.6+8.4) 230 (0+734) 
Barley 227 (173+311) 120.1 (0+371) 


4. CONCLUSIONS 


At development of a methodology to assess critical or no-effect doses for biota 
as a whole or for particular groups of flora and fauna, for example, 
domesticated plants in this work or vertebrates, invertebrates and plants in 
[14], there are many questions to be answered. It is commonly agreed that 
protection should be provided not for individuals but for populations and species 
to preserve biodiversity in ecosystems. Trying to reveal biological effect in 
species most suitable for verification of radiation safety standards, we have 
found that traits of productivity and reproduction success in cultivated plants 
appear to be more sensitive to radiation than survival endpoint. An adaptation 
of the ERA methodology for the environment radiation protection [5, 6, 14] is a 
promising idea in our opinion, but some of its key techniques need a further 
testing. Thus, it is important to find an effective way of data treatment to obtain 
reliable estimates. The present study gives preliminary findings contributing to 
development of criteria to assure quality of input data and justify a dose-effect 
dependence shape for the environmentally relevant effects. In particular, it 
shows that an absence of practice to apply standard ecotoxicity tests in 
radiation toxicology leads to an inadmissible waste of experimental data (80% 
and more) when the ERA approach is followed. At acute exposures, the 
proposed approach works well, and critical doses obtained from different 
models are comparable. In the most important from the viewpoint of the 
environment protection case of chronic exposure, however, conclusions made 
from different approaches to data treatment can differ significantly. At this, 
critical dose rates estimated from the linear model seem to be less scattered 
and better defined than at the logistic model utilization. However, a lack of 
proper data makes any differences insignificant. It is obvious that a further data 
gathering as well as development and testing of techniques for their treatment 
to get justified estimates of critical doses and dose rates for biota are needed. 
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